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Abstract 

We provide an elementary proof to a conjecture by Robinson that multiples of (powers of) 
primes in the Somos-4 sequence are equally spaced. We also show, almost as a corollary, for 
the generalised Somos-4 sequence defined by r n +2T„_2 = ar„+iT„_i + /3t% and initial values 
ti = T 2 = T3 = T4 = 1, that the polynomial T n (a, ft) is a divisor of T n+ k(2n-s)(ct, P) for all 
n,k £ Z and establish a similar result for the generalized Somos-5 sequence. 


1 Introduction 

An arithmetic divisibility sequence is neither a divisibility sequence, nor an arithmetic sequence. 
We define it to be a sequence {f n }ne z such that there is a function d : Z i-> N, called the common 
difference function, such that 


Vn G Z : d{n) | (to - n) => /„ | f m . 

Clearly, if for some sequence {f n }nei prime powers appear in arithmetic progressions that sequence 
is an arithmetic divisibility sequence. Furthermore, we then have the following nice property: for 
all to = n + kd(n) and r = m + sd(m), with n, k, s £ Z, there is an l £ Z such that r = n + ld(n). 
Primes do not have to appear in arithmetic progressions, e.g. consider the periodic arithmetic 
divisibility sequence 

12 if n = 0 mod 6, 

1 if n = 1, 5 mod 6, 

6 if n = 2,4 mod 6, 

4 if n = 3 mod 6. 

where we have even made sure that f n \ f m => d(n) | (to — n). The above mentioned nice 
property also holds if the common difference function is polynomial. This follows from the fact 
that p{x) divides p(x + p(x)q(x)) for all polynomials p(x), q(x), which can be seen by substituting 
y = x+ p[x)q(x) in y-x \ p{y) - p(x). 

Sequences {r n } defined by 

Ai+fc—2 Ai—2 = O^'Xn+k —3 Ai—1 d“ /?V n -f_/ c _4 T n (1) 

possess the Laurent property, that is, all their terms are Laurent polynomials in their initial 
values. The equations m arise as a special case of the Hirota-Miwa equation, whose integrability 
condition is equivalent to Laurentness [2D] - For k = 4 and k = 5 these equations are called 
Somos-4 and Somos-5. The Laurent property implies in particular that if the initial values are 
1, and a,j3 are integers, both recurrences define integer sequences. Considering a and /? to be 
variables, we obtain polynomial sequences. The Laurent property was introduced by Hickerson, 
cf. [23] . to prove the integrality of the next sequence, Somos-6, which does not have the form m 
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but is a special case of the discrete BKP equation [23 EU. Initially the integrality, and secondly 
the Laurent property itself, seemed to be curious properties m- By now there is extended 
literature on equations that possess the Laurent property, have Laurentness, or display the Laurent 
phenomenon, n e i m eh m m uni m , and the property is well understood in terms of cluster 
algebras [_?j, which have had a profound impact in diverse areas of mathematics, cf. [ 6 |. The Somos- 
k = 4, 5 sequences also satisfy a strong Laurent property, e.g. for Somos-4 the terms are Laurent 
in n, polynomial in T 2 ,..., t*,, and polynomial in X, an additional (invariant) rational function of 
the initial values [15] . Therefore, integral sequences may be obtained for other than unit initial 
values. The Somos-4 sequence with a = — 1, (3 = 2 and initial values t\ = 72 = 1, T 3 = 2, T 4 = 3 
is an integer sequence which extends the Fibonacci numbers; we have r n = (—l) n+ 1 r_ n and 
= 7n— 1 4“ ^~n— 2- 

Few results on divisibility properties of the Somos-(fc = 4, 5) sequences (with a = f3 = 1 and 
7 ~i = • • • = Tfc = 1) are known. In [3] it is shown that every term beyond the fourth of the 
Somos-4 sequence has a primitive divisor, i.e. a prime which does not divide any preceding term. 
Jones and Rouse m show that the density of primes dividing at least one term of the Somos-4 
sequence is 11/21. In |18 the authors establish that the terms of Somos-4 (with a = (3 = 1) 
are irreducible Laurent polynomials in their initial values and pairwise co-prime. Robinson [23] 
showed that, for k = 4,5, the 7-th and j-th terms of the Somos-fc sequence are relatively prime 
whenever |i — j\ < k. He infers that both sequences are periodic modulo m for every m £ N. In 
[23 section 4] he presents results about Somos-4 and Somos-5 obtained by calculation, for which 
he did not have general proofs. One observation he made is that if a prime p divides any term of 
Somos-fc, then the multiples of p are equally spaced. In proof he added that Clifford S. Gardner 
proved this result, but a reference was not provided. This conjecture and generalisations have been 
proven by C.S. Swart in her thesis [26] . She also proved, for more general Somos-4 sequences, that 
either all multiples of a prime p are divisible by exactly the same power of p , or 3r £ NVA; >rfM 
some term is divisible by p k exactly, cf. section [ 8 ] The following problem is still open: to prove 
that for particular Somos sequences for all but finitely many primes p the multiples of p are not 
divisible by exactly the same power of p. 

In this short paper we provide a simple proof, in section [2 that for all primes pgP the set 
{n £ Z : p k | r„} has either less than 2 elements, or its elements form a complete arithmetic 
sequence. This also holds for rational Somos-fc = 4, 5 sequences, for primes that do not appear in 
the denominator of the invariant function, as long as the initial conditions are pairwise co-prime. 
It also holds true in elliptic divisibility sequences, see section [3] We employ an elementary result 
from number theory, which we haven’t found in the literature, namely that the terms of a subset 
S C Z with the property that 2s — t £ S for all s,t £ S, form an arithmetic sequence. A proof 
of this result is presented in section [2] Combining the result on the appearance of primes with a 
symmetry of the equation and taking unit initial values we find Somos-fc = 4, 5 sequences which 
are arithmetic divisibility sequences with common difference function dfc(n) = 2 n — k — 1 , i.e. we 
prove the polynomial divisibility 

T n (a , /3) | T n+ i dk{n) (a, (3) Vn, l £ Z. 

We briefly investigate different initial values which give rise to Laurent sequences with the same 
divisibility property, as well as initial values that give rise to polynomial sequences with slightly 
different divisibility properties. We note that if one could prove for odd primes p that the multiples 
of p in the Somos-4 & 5 sequences are not divisible by exactly the same power of p then we would 
have the following. 

Conjecture 1. Let r n be the terms of the Somos-k sequence with k = 4,5, a. = /3 = 1, ti = • • • = 
Tfc = 1, d = 2n — k — 1 and q = r n when (k + 1) \ n or q = T n /2 when (k + 1) | n. Then 

q m+1 | n^l = n+ ( i -^— + kq m )d. 

A similar result for the extended Fibonacci sequence would provide an alternative proof for 
the fact that the (m + l)-st power of the ?r-th Fibonnaci number divides the (n/™)-th Fibonnaci 


2 








number, as was shown by M. Cavachi 3], at the age of 16. We note the case n = 3 is exceptional; 
here 2 m + 2 | / 2 » 3 . 


2 Sets of differences 

In [2] a set of differences is defined to be a subset S of N such that s,t £ S, s > t => s — t £ S. 
A useful result is that the elements in sets of differences are multiples of the smallest element. 
Extending the definition to subsets of Z a similar statement may be used to provide a quick proof 
that elliptic divisibility sequences are divisibility sequences. 

For our purposes we define a subset S of Z to be a modified set of differences if s, t £ S => 
2 s-t £ S. 

Lemma 2. A modified set of differences has less than two elements, or its terms form a complete 
arithmetic sequence. 

Proof. Let s, t be elements of a modified set of differences V. We first show, by induction, that 

(k + l)s — kt £ V, and (k + 1 )t — ks £ V, (2) 

for all k £ Z. By definition it is true for k = 1. Assuming the statement we obtain that 
2 s — ((k + l)s — kt) = kt — (k — l)s and 2 1 — (( k + l)s — kt) = (k + 2 )t — (k + l)s are in V, as well 
as the same expressions with s, t interchanged. 

If s = t all the above elements coincide, the set may only have one element. If s and t are 
distinct we may take s < t in which case ( k + 1 )t — ks with k £ N provides plenty of positive 
elements. Let a,b > a be the smallest and the next smallest positive elements of V and denote 
d = b — a. As a — d = 2a — b £ V and a is the smallest positive element we have a < d. Taking 
s = a and t = 2a — b in the first element of (0) we find that kd + a £ V for all k £ Z. We show 
the converse, s£f => d \ (s — a). 

For any s £ V we may write s = qd + r with 0 < r < d. If r = a then we are done. Let 
t = qd + a £V. Assume r < a. The element z = (k + l)s — kt = s — k(a — r), with k = L^zyJ - is 
in V and 0 < z < a — r which is impossible, as a is the smallest positive element of V. Assuming 
r > a, we have, with k = J > that p = (k + 1 )t — ks = t — k(r — a) is in V and 0 < p < r — a. 
We also have, replacing k with k — 1 , that q = p + (r — a) £ V. As p < d < b we need p = a. But 
with p = a we have a<q = r<d<b which contradicts that b is the second smallest positive 
element. □ 


3 Elliptic divisibility sequences, and integrality 

Elliptic divisibility sequences (EDS) were introduced by Morgan Ward \T7\ [55] as sequences of 
integers {a n }ffL 0 , that satisfy, for all m > n > 1 , 


Q j m-\-nQ j m—n 


^n+l^m Q j nQ j m -\-1 


(3) 


and n \ m =£> a n | a m . He calls a sequence proper if do = 0, cq = 1 , a\ + a 3 0, and shows that 
a proper solution to ( 0 ) is an EDS if and only if 02,03,04 are integers and a 2 \ 04. Ward first 
shows by induction, that all terms are integers, and then by another induction step, the divisibility 
property. 

Ward’s proof of integrality is generalised by Hone and Swart m, who proved the following 
strong Laurent property for Somos-4: the terms r„ are polynomials in a, f3, t^ 1 , t 2i 73, T4, and Z, 
where 


1 = a 2 + f}T, 


t 2 r| + a(r| T 4 + rir|) + /3 t|t| 


(4) 


TlT 2 7-3 t 4 
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Therefore, if n = ±1 and I is integer, the sequence consist of integers. We remark that taking 
n = 2 in © gives us a special Somos-4 sequence, 

Om+ 2 Qm -2 = — 1 (5) 

Taking a = a|,/3 = — 0103 , + = a ,}^ =1 and of = 1 in d4]) we find I = —< 34/02 , whose integrality 
implies that the sequence {a m } consists of integers. 

For Somos-5 we have the following strong Laurent property [T5J Theorem 3.7]: The terms r n>0 
of Somos-5 are polynomial in a, /?, iq, + 1 , 73 , T 4 , T 5 , and J , where 

J = f3 + aS S = ( T i r 5 + orf )( r i T l + r 2 +) + + 2 T 3 T 4 

T 1 T 2 T 3 T 4 T 5 

The sequence is an integral sequence if ri, 12 G {±1}, I is integer and one can find three consecutive 
integers preceding two units. 

An EDS uniquely extends to a sequence over Z. Taking m = 2 in (0 one finds ao = 0. Taking 
m = 1 in (0 we find 030-1 = — 03030 ^ => a_i = —1 (assuming 03 7 ^ 0). Taking m = — 1 and 
n = k — 1 in @ gives ak- 2 d-k = —Ofc_ 20 fca _ 1 => a_fc = —Ofc (assuming Ofc -2 7 ^ 0). For k = 2 this 
doesn’t work, but there exist a determining equation, e.g. take m = 0 in (f5l) . One way to deal 
with the occurrence of zeros other than ao is to take the initial values as parameters (let 04 be 
a multiple of 02 ), generate a polynomial sequence, and then specialise. See 15, Appendix A] for 
another discussion on zeros. 

Taking the initial values as parameters (or assuming them to be co-prime) one can inductively 
show that gcd(a„, a n +i)=l for all n > 1. The divisibility property follows by showing that 
Vk = {to G Z : a,k | a m }, with k > 1, is a set of differences with at least two elements. To show 
this we employ a second family of recurrences 


CL 1 Ci2 dm-\-n -\-1 —n 


Cl n CLm — l dri—l &rn 

^ 71 + 2 ^ 771+1 ^ 71 + 1 ^ 771+2 


( 6 ) 


Taking m = s and n = t — s in @ gives 

2 2 
0t02s-t = 1t-s a s-l a s+l — at-s-iat-s+lCi s - 

If s G Vk then s ± 1 ^ 14, and thus Ofc | af_ s . Taking m = s and n = t — s — lin © gives 


CL\CL2CltCL2s —£+1 — CLt—s— l^t—l^s+2 &t—s—2^t— s+l^s^s+1 1 

and so Ofc | at- s -iat- s ■ Together with cit- s -± and dt- s being co-prime we find that s,t G T4 => 
s — t G Vk- As both 0 and k are elements of 14 we have 14 = kZ. Thus, the polynomial sequence 
{a n } is an arithmetic divisibility sequence with common difference function 


d(n) 


n n + 0, 
oo n = 0. 


For special initial values the value of d at 0 can be finite. 

Ward does not give much detail about the consistency of the family of recurrences [3] He does 
however provide an explicit solution for a n in terms of the Weierstrass sigma function, and both 
© and © are direct consequences of the corresponding three-term relation, cf. M- For an 
algebraic approach we refer the reader to [2Tj. 

The fact that an EDS {a ra } is a divisibility sequence does not imply that multiples of powers 
of primes are equally spaced. However, for an EDS whose initial values satisfy gcd(a n , a n +i)=l 
the above argument shows that V = {n G Z : p k \ a ra }, with p G P, k G N is a modified set of 
differences with at least two elements, and this implies the following theorem. 

Theorem 3. For an EDS in which subsequent (initial) values are co-prime the multiples of powers 
of primes are equally spaced. 
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4 Compagnion elliptic divisibility sequences 


Taking unit initial values for Somos-4, 


n = r 2 = r 3 = r 4 = 1, (7) 

we have T = 1 + 2a + /3, and X = a 2 + /3T = (a + /3 ) 2 + /3. We define JT5J Definition 1] an EDS 
by ([5]) and the initial values 

a 2 = —y/a, 03 = —/3, 04 = \/aX ( 8 ) 

so {a n } satisfies the same recurrence as {r„}. The sequence {a n } is the companion EDS for Somos 
4, that is, the following families of recurrences are satisfied m corollaries 1.2, 1.3], 




Q"n7~m —1 


and 


tti(X2T m -|_ n -|-iT m _ n — 


U j n'Tm—l Q'n—l'Tm 
Q j n+2'T~m+l ^n+l^m+2 

Proofs of these facts can be found in [21 UQ > cf. [13] . 

To describe the companion EDS for Somos-5 with initial values 


Tl = T2 = T3 = T4 = T5 = 1, 


(9) 

( 10 ) 


( 11 ) 


we introduce an alternating sequence of functions 


hi 


2a + /3 l = 0 mod 2, 
a + 1 1 = 1 mod 2 , 


whose product equals the above mentioned invariant J = hihi + 1 , see [151 Proof of Theorem 3.7]. 
The companion EDS is then given by 

«i = 1 , a 2 = y/hm+m a 3 = a, a 4 = -/3a 2 , 


and 

Q'k-\-2Q'k—2 = 

We note that these are actually two companion sequences, one for each value of n + m mod 2, and 
we have 

a 2 k {m + n = i) _ hi 
a 2 k (m + n = j) hj’ 

so any polynomial divisor of afc(O) is also a divisor of afe(l) and visa versa. The terms of the 
Somos-5 sequence satisfy the same families of recurrences mm , see also m corollary 2 . 12 ] and 
EH Section 7]. 


5 Relative primeness 

Both a and /3 are not divisors of r ra . This can be seen by taking a = 0 or ft = 0. The corresponding 
solutions, for Somos-4, are T n = f3 kn and r n = a ln , where k n and l n satisfy the linear recurrences 

kn- f -2 = 2 k n k n —2 T 1, ln+2 = ln+1 T In— 1 In—2 T 1- 

In particular, they do not vanish. If a and /3 have a divisor in common its multiplicity will grow as 
f [24l A249020]). Considering the terms r n as polynomials in a and /3, we can adjust the argument 
of Bergman, cf. m, to prove that any four consecutive terms of our polynomial Somos-4 sequence 
are pairwise co-prime. The initial terms 0 are co-prime. Assume that {T n +i }} = _ 2 are pairwise 
co-prime, and let p be an irreducible divisor of r n+2 with positive degree. Obviously p does not 
divide a or (3, and therefore p | a n if and only if p \ a n -1 or p | a n + 1 . By hypothesis this does not 
happen. Similarly for Somos-5 with pairwise co-prime initial values any five consecutive terms are 
pairwise co-prime. 
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6 Divisibility in Somos-4 and Somos-5 

Theorem 4. The Somos-k sequences, for k = 4, 5, with initial values T\ = • • ■ = = 1 are 

arithmetic divisibility sequences with common difference function d{n) = 2n — k — 1 . 

Proof. We first assume that the initial values are co-prime and show that the set W q = {n € Z : q \ 

r„} is a modified set of differences. Here r n is a solution of either Somos-4 or Somos-5, and we will 

take q to be any polynomial in the ring Z[a,/ 3 ,I,rj ± 1 ,T 2 ,T 3 ,T 4 ], or Z[a, j3, J, 73 , 74 , 75 ], 

respectively. Of course, one might take q = Tk, or let q = p h where p does not occur in a 

denominator. Taking m = s and n = t — s in J9| gives us 

2 2 
TtT2s—t = Oj^_ s Ts— lTs-|-i CLt—s—l^t— s+l'T’g • 

If s, t £ then s ± 1 ^ W q and hence q \ a 2 _ s . Taking m = s and n = t — s — lin m yields 
q | at- s at- s - 1 , and hence q \ at~ s ■ Taking m = s and n = t — s in Hoi gives us 

aia 2 Tt+iT 2 s-t = a t _ s a t _ s+ iT s _iT s+2 — at- a -iat- s + 2 T s T s +i. 

As q j aia 2 T t+1 we obtain 2s — t £ W q which is what we wanted to show. By lemma [2] it follows 
that W q has less than 2 elements or there exist n, d £ N such that m £ W p => d \ m — n. 

We now specialise to initial values ri = • ■ ■ = Tfc = 1. Due to the symmetry of the recurrences, 
T n +i -O- r n -i for all l, we obtain T— n+ (k+i )/2 = T n +tk+ 1)/2 for all n. In terms of d = 2n — k — 1 we 
have r n -d = r n . As the degree of rfe, with 0 < k < n, is smaller than the degree of r„, we have 
r n \ Tk and hence d is the common difference. □ 


7 Equivalent sequences 


It is possible to choose other initial values with the same symmetry. For Somos-4, taking t[ = 
—t 2 = — 7-3 = 74 = 1 one gets the polynomial sequence T' n {ot,f$) = (—l)^ n ^ 2 W n (—a,/3), where 
T n (a,j3) is the sequence obtained from Ti = r 2 = 73 = 74 = 1. More generally, starting from 
t[ = 7-4 = 1 and t 2 = t 3 = 7 we find a sequence in Z[a, 7 ±1 ], 


T~n(WP) 


T- ra (7 3 a,7 4 /3) 
-y(«-1)(«-4)/2 ' 


( 12 ) 


If ft = 7 2 the sequence is polynomial, and for n > 4, with d = 2n — 5, 


n = 1 mod 3 : T' n \ T' m 4=> d \ (m — n), 
nj^l mod 3 : r' n \ T' m 4=> d | (to — n) and 


^ 1 mod 3. 


Furthermore, starting from r( = 74 = d and T 2 = T 3 = 7 we find a sequence in Z[a, /3, 7 ±:l , d^ 1 ] 

£(n—2)(n—3)/2 / /n/ \3 /n/\4 




(n—l)(n—4)/2 




Both (J 1 2|) and m have the same divisibility properties as T n (a,j3). 
For Somos-5, starting from t[, ... ,tI = a, b, c, b, a we have 


TnfaP) = 




/ 

( c \ 2 

/ c. \ 3 \ 

-T„ 

( £ ) a, 

(- p) 

\ 

\aJ 

\aJ I 


(13) 


(14) 


where 

4 n 2 3n 17 — (—1)" D 1 + (-1)" ^ n 2 3n 13 + 3(-l) n 

n ~T _ T + 8 ’ 2 ’ n ~T _ T + 8 


We get a polynomial sequence from initial values 74 = 75 = 1, 12 = T 4 = b, and t 3 = a. Here we 
find the divisibility, with d — 2n — 6, r n | r m 47 d | (to — n) for n > 5, r„ | r m Sn|m for n = 2,3. 

We note the sequences r and t' in this section are equivalent sequences in the sense of (25j, 
Sect. 6.3], and that the slightly different divisibility properties are due to initial values having a 
common factor. 
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8 Robinson’s observations 


We state here some observations made in [23] and their current status. Robinson observed that 

1 . the multiples of primes are equally spaced [such a prime is called regular ] 

2 . the gap (the common difference) is never much larger than p 

3. if p occurs then p 2 occurs, and its gap is p times the gap of p 

4. if p 1 is the smallest occurring power, the gap of p l+l will be p l times the gap of p l 

Much of this is now understood, but not all. Hone and Swart hi uncoil (as well as Naom Elkies, 
David Speyer, and Nelson Stephens in unpublished work, cf. [22]) have shown that the terms r„ 
of a Somos-4 sequence correspond to rational points Q + [n\P on an associated elliptic curve E. 
Christine Swart, in her thesis [26] . studies elliptic curves over Z p r- due to an equivalence 


T n = 0 mod p r -£=> Q + [n]P = O p r , 


if p r occurs, the gap of p r equals the order N r of the (non-singular) point P in E(Z p r). Swart has 
proved that either all powers p k are regular, or all multiples of p are divisible by exactly the same 
power of p [2B] Thm 7.6.6]. She obtained the structure of the gap-function [26] Thm 7.6.7] and 
explained (and improved) Robinson’s bound on the gap exploiting the Hasse bound which bounds 
a multiple of the order of P, i.e. the number of points on F(F p ) within 2 ^Jp of p + 1 [26] Thm 
7.6.5]. Armed with the above mentioned theorem, the gap can be calculated from the order of P 
on E(F P ). 

We conclude with a couple of examples of interest, using curves E and points P as given in 

m- 

• Somos-4 with a = /3 = t± = • • • = = 1. The prime 2 divides r m if and only if 5 divides 

to. Higher powers of 2 do not appear (the sequence mod 4 is periodic with period 10 and 
does not contain 0). It seems that for all p / 2 such that 3n g N : p | r n , all powers 
p k are regular, and N r = p r ~ 1 Ni, where Ni is the order of the point (1,1) on the curve 
y 2 = 4(x 3 — x) + 1 modulo p. To prove conjecture Q] it suffices to show that all powers of 
odd primes occurring are regular, one does not need w = 1 for all p. 


• Somos-4 with a = —1, (3 = 2 , (n,--- , T 4 ) = (1,1,2,3) (which extends the Fibonacci 
sequence). For p = 2 we have w = 1 and v = 2, for p = 3 we have w = 1 and Ni = 4. It 
seems that all powers of occurring primes are regular, N r = p r ~ w N\ where N\ is the order 
of the point yj— 1) on the singular curve 

2 (6a; — 5)(5 + 12a;) 2 

V ~ 216 


over the field F p [-\/—1]. Again, we have not found a value of w different than 1 . 

• Somos-4 with a = 4, (3 = 9, T\ = T 4 = 1, r 2 = r 3 = 3. We have 3 k \ r m k = 1, 3 f to. 

Taking p = 5 we find N r >3 = 5 r_3 Wi where N\ =7 is the order of the point (55750/243,2) 
on the curve given by 


y 2 = 4ar - 


12428112196 
19683 3 


1385503884676628 

14348907 


over F 5 . 


• Somos-4 with a = 2, (3 = 5 , (ti, . .., 74) = ( 1 , 3, 2, 5 ). We have 7 k \ T m <t=> k = 2, Ni \ m, 
where Ni = 10 is the order of P = (223081/21600, \/2) on 

2 „ , 48492460561 10678311547192441 

1 / = 4:j: —- x -b- 

y 38880000 1259712000000 

over F 7 or over Ft[\/ 2 ]. 
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